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Abstract 

In this paper we derive the Clauser-Horne (CH) inequality for the full electron counting statistics 
in a mesoscopic multiterminal conductor and we discuss its properties. We first consider the 
idealized situation in which a flux of entangled electrons is generated by an entangler. Given a 
certain average number of incoming entangled electrons, the CH inequality can be evaluated for 
different numbers of transmitted particles. Strong violations occur when the number of transmitted 
charges on the two terminals is the same (Q\ = Q2), whereas no violation is found for Qi ^ Q2- 
We then consider two actual setups that can be realized experimentally. The first one consists of 
a three terminal normal beam splitter and the second one of a hybrid superconducting structure. 
Interestingly, we find that the CH inequality is violated for the three terminal normal device. 
The maximum violation scales as 1/M and 1/M 2 for the entangler and normal beam splitter, 
respectively, 2M being the average number of injected electrons. As expected, we find full violation 
of the CH inequality in the case of the superconducting system. 
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I. INTRODUCTION 



Entanglement [l| denotes the nonlocal correlations that exist, even in the absence of direct 
interaction, between two (spatially separated) parts of a given quantum system. Since the 
early days of quantum mechanics, understanding the phenomenon of entanglement has been 
central to the understanding of the foundations of quantum theory. Besides its fundamental 
importance, a great deal of interest has been brought forth by its role in quantum informa- 
tion Entanglement is believed to be the main ingredient of computational speed-up in 
quantum information protocols. 

Most of the work on entanglement has been performed in optical systems with photons 0] , 
cavity QED systems Q and ion traps Only recently attention has been devoted to the 
manipulation of entangled states in a solid state environment. This interest, originally moti- 
vated by the idea to realize a solid state quantum computer 0, 0] , has been rapidly growing 
and by now several works discuss how to generate, manipulate and detect entangled states in 
solid state systems. It is probably worth to emphasize already at this point that, differently 
from the situation encountered in quantum optics, in solid state system entanglement is 
rather common. What is not trivial is its control and detection (especially if the interaction 
between the different subsystems forming the entangled state is switched off). 

Despite the large body of knowledge developed in the study of optical systems, new 
strategies have to be designed to reveal the signatures of non-local correlations in the case 
of electronic states. For mesoscopic conductors, the prototype scheme was discussed in 
Ref. In this work it has been shown that the presence of spatially separated pairs of 
entangled electrons, created by some entangler, can be revealed by using a beam splitter 
and by measuring the correlations of the current fluctuations in the leads. Provided that 
the electrons injected are in an entangled state bunching and anti-bunching behavior for 
the cross-correlations of current fluctuations are found depending on whether the state is a 
spin singlet or a spin triplet. Not only the noise, but the full counting statistics is sensitive 
to the presence of entanglement in the incoming beam The distribution of transmitted 
electrons is binomial and symmetric with respect to the average number of transmitted 
charges. Moreover, this is important for the problem studied in the present work, the joint 
probability for counting electrons at different leads unambiguously characterizes the state 
of the incident electrons if one uses spin- sensitive electron counters. In this case the joint 
probability cannot be expressed as a product of single-terminal probabilities. 

Given the general setup to detect entanglement an important issue is to understand how 
to generate it. This has been discussed in several papers. Most of the existing proposals 
are based on the generation of Bell states by means of electron-electron interaction. This 
can be achieved through superconducting correlations I(| in hybrid normal - superconduct- 



ing [111. |12|. Il3l. |14|] and superconductor - carbon nanotubes systems [15l. |16|. quantum dots 



in the Coulomb blockade regime [17[ or Kondo-like impurities [18|. Then, by using energy 
or spin filters, the two electrons forming the Bell state are separated. The entanglement 



can be created in the spin or in the orbital [14| degrees of freedom. Very recently, as it 



is also discussed in Section IIIIB| it was shown that in a mesoscopic multi-terminal con- 
ductor entanglement can be produced also in absence of electron interaction [l9|. Besides 
electrons, it ispossible to produce entangled states with Cooper pairs in superconducting 



nano circuit s |20| or by coupling a mesoscopic Josephson junctions with superconducting 



resonators [21|, |22j, |23L|24 



Since Bell's work [251 ] . it is known that a classical theory formulated in terms of a hidden 
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variables satisfying reasonable condition of locality, yields predictions which are different 
from those of quantum mechanics. These predictions were casted into the form of inequal- 
ities which any realistic local theory must obey. Bell inequalities have been formulated for 
mesoscopic multi-terminal conductors in Refs. [IH, 26|, 27j in terms of electrical noise corre- 
lations at different terminals [28[ . A test of quantum mechanics through Bell inequalities in 
mesoscopic physics is very challenging and most probably it would be rather difficult, if not 
impossible, to get around all possible loopholes. Although solid state systems are not the 
natural arena where to test the foundations of quantum mechanics, it is nevertheless very 
interesting to have access, manipulate and quantify these non-local correlations. 

In this work we derive a Bell inequality for the full electron counting statistics and discuss 
its properties. The formulation we follow is based on what is known as the Clauser-Horne 



(CH) inequality [29|, [30|. We shall show that the joint probabilities for a given number of 
electrons to pass through a mesoscopic conductor (in a given time) should satisfy, for a 
classical local theory, an inequality. 

The paper is organized as follows: in the next Section we motivate our approach to the 
problem, derive the CH inequality and express the joint probabilities needed in the CH 
inequality in terms of the scattering properties of the mesoscopic conductor. Section II III 
is devoted to the discussion of the results. We first consider the idealized situation where 
an incoming flux of fully entangled electrons is injected into the mesoscopic region. Then 
we move on to analyze actual setups. Interacting electrons are not necessary to have an 
entangled state, we show that a three terminal normal device is enough to lead to violation 
of the CH inequality. For completeness we also consider the case where entanglement is 
produced by Andreev reflection. In the last Section we present the conclusions and a brief 
summary of this work. 



II. CH INEQUALITY FOR THE FULL COUNTING STATISTICS 

Electron Full Counting Statistics (FCS) refers to the probability that a given number of 
electrons has traversed, in a time t, a mesoscopic conductor. In the long time limit the first 
and the second moment of the probability distribution are related to the average current 
and noise, respectively. The reason for which we resort to FCS for analyzing electronic 
entanglement in a solid state environment resides in the fact that electrons in a conductor 
are not necessarily sufficiently separated from one another for coincidence counting to make 
sense, like in optical systems. Furthermore, the measurement of single coincidence events 
in electronic solid state systems does not seem realizable at present. Zero-frequency noise 
accounts for long time correlations and we do not expect it to be in general sensitive to 



coincidence measurements (see however the discussion in Ref. [14J for the limit of small 
transmission rates). From these premises we suggest that FCS is a natural candidate to 
formulate a Bell-type inequality for electrons in mesoscopic conductors. In the case where 
only two entangled electrons are injected, we find a situation similar to that with photons. 
More generally we discuss the case where a large number of electrons have been injected. 



In its original version [25[, the Bell inequality was derived for dicotomic variables. Here 



we consider the more general formulation due to Clauser and Home [29|]. We consider the 
idealized setup, illustrated in Fig{TJ which consists of the following parts. On the left we 
place an entangler that produces 2M electrons in a spin entangled state (in Section II I II two 
different situations for the implementation of the entangler are discussed). Two conductors, 
characterized by some scattering matrix, connect the terminals 3 and 4 of the entangler 



3 



with the exit leads 1 and 2 so to carry the two particles belonging to each pair into two 
different spatially separated reservoirs. The electron counting is performed in leads 1 and 
2 for electrons with spin aligned along the local spin-quantization axis at angles 9\ and 82- 
Detection is realized by means of spin-selective counters, i.e. by counting electrons with 
the projection of the spin along a given local quantization-axis. In analogy with the optical 
case we say that the analyzer is not present when the electron counting is spin-insensitive 
(electrons are counted irrespective of their spin direction). 

In Section THAI we present the derivation of the CH inequality for the FCS and in Section 
III Bl we resume, for completeness, the relation between FCS and the scattering matrix S. 



A. Derivation of the CH inequality 

The basic object for the formulation of the CH inequality is the joint probability P(Qi, Q2) 
for transferring a number of Qi and Q2 electronic charges into leads 1 and 2 over an obser- 



vation time t. We follow closely the derivation given in Ref. [30|. Our starting point is the 
following algebraic inequality 

— 1 < xy — xy' + x'y + x'y' — x' — y < (1) 

which holds for any variable < x,y,x',y' < 1. Let us now introduce explicitly a set of 
hidden variables r which take values in a space T. We assume that the incoming entangled 
electron states are described by r in all the details necessary to determine the probability 
distributions P(Q a , r) for transferring a number of Q a electronic charges into lead a = 1, 2. 
By imposing that the hidden variable theory is local, it follows that the joint probability 
can be expressed in the following form: 

P(Qi,Q 2 ) = J M(r)P(Q u r)P(Q2,T)dr, (2) 

where Ai(r)dr defines a probability measure on the space T. The physical meaning of Eq.([2]) 
is straightforward: it states that the probability distribution on lead a does not depend on 
the probability distribution on the lead /?. 

We now introduce P 6l,92 (Qi, Q2) as the joint probability for transferring Q 1 and Q2 
electronic charges when both analyzers are present, while P ~(Qi, Q2) an d P~' e ' 2 {Qi,Q2) 
are the corresponding joint probabilities when one of the two analyzers is removed. If the 
condition 

P 6a (Q a ,r)<P(Q a ,r) (3) 

(known as no- enhancement assumption ) is verified, it is possible to identify the variables 
appearing in Eq.([T]) as follows: 



P dl (Qi,r) P* 2 (Q 2 ,t 

x = — r V 



X 



P(Q u r) » P(Q*,t) 
P e KQi,r) ,_ P e ^Q2,T) 
P{Qx,t) V ~ P(Q 2 ,t) 



(4) 



P e (Qa,i~) being the single terminal probability distribution in the presence of a analyzer. 
Eq. ([1]) can then be rewritten in terms of probabilities by multiplying each side of the equa- 
tion by P(Qi, t)P(Q2, r)A4(r)dr and integrating over the space T. Finally the following 
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inequality is obtained 

Sen = P eifi KQi,Q2)-P^HQi,Q2) + P^ e KQuQ2) + P e ^KQuQ2) 

- P e 'i'-(Qi,Q2)-P~> d2 (Qi,Q2)<0. (5) 

Eq.® is the CH inequality for the full counting statistics 46[, holding for all values of Qi 
and Q2 which satisfy the no-enhancement assumption. We stress that the no-enhancement 
assumption, upon which Eq. (jSJ) is based, it is not satisfied in general like its optical version. 
The quantities that we have to compare are probability distributions, so that Eq. (jSJ) must be 
checked over the whole range of Q. For a fixed time t and a given mesoscopic system, hence 
for a given scattering matrix and incident particle state, the no-enhancement assumption is 
valid only in some range of values of Q. In particular, different sets of system parameters 
correspond to different such ranges. The quantity Sch in Eq.® depends on Q\ and Q2 so 
that the possible violation, or the extent of it, also depends on Qi and Q2- Given a certain 
average number M of entangled pairs that have being injected in the time t, one can look 
for the maximum violation as a function of the transmitted charges Qi and Q2- 



B. Scattering approach to the full counting statistics 



The joint probabilities appearing in Eq. (l5"]) can be determined once the scattering matrix 
5" of the mesoscopic conductor is known. The FCS in electronic systems was first introduced 



by Levitov et al. in Ref. |32l . |33| in the context of the scattering th eory and later on the 
Keldysh Green function method [HJ to FCS was developed in Refs. 35] (for a review see 
Refs. j36[). In this paragraph we briefly describe how the FCS is formulated for a mesoscopic 
conductor in the scattering approach. Within this framework, the transport properties of 
a metallic phase-coherent structure attached to n reservoirs are determined by the matrix 
S of scattering amplitudes [37j|. Such amplitudes are defined through the scattering states 
describing particles propagating through the leads. For one dimensional conductors, for 
example, the scattering state arising from a unitary flux of particles at energy E originating 
in the i-th reservoir reads 



ipi{x) 



n{E)e 



-iki(E)x 



for the i-th lead, and 



tji(E)e- 



-ikj(E)x 



(6) 



(7) 



for the j-th lead, with j 7^ %. Here r^E) is the reflection amplitude for particles at energy E, 
wave vector ki(E) and group velocity ih(E) and tji(E) is the transmission amplitude from 
lead % to lead j. Note that |rj| 2 is the probability for a particle to reflect back into the i-th 
lead and \tji\ 2 is the probability for the transmission of a particle from lead i to lead j. In 
the second quantization formalism, the field operator ipj a (x,t) for spin a particles in lead j 
is built from scattering states and it is defined as 38 



dE 



a ja (E)e 



,a{E)e 



(8) 
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where aj a {E) (4>j a (E)) is the destruction operator for incoming (outgoing) particles at energy 
E with spin a in lead j. These operators are linked by the equation 



/ 01T \ 



,5' 



/ «1T \ 

ail 

\ i / 



(9) 



V = / 

and obey anti-commutation relations 



(10) 



In the case of two and three dimensional leads one can separate longitudinal and transverse 
particle motion. Since the transverse motion is quantized, the wave function relative to 
the plane perpendicular to the direction of transport is characterized by a set of quantum 
numbers which identifies the channels of the lead. Such channels are referred to as open when 
the corresponding longitudinal wave vectors are real, since they correspond to propagating 
modes. Note that the case of a single open channel corresponds to a one dimensional lead. 

Let us now turn the attention to the probability distribution for the transfer of charges. 
Following Ref. 3^], within the scattering approach the characteristic function of the proba- 
bility distribution for the transfer of particles in a structure attached to n leads at a given 
energy E can be written as 



XE(\ h \.l) = (] 



e x ji 



j=l,n 



n 

j=l,n 



it 

° e 



-initio 



(11) 



where the brackets (...) stand for the quantum statistical average over the thermal distri- 
butions in the leads. Assuming a single channel per lead, Nj? -> is the number operator for 

incoming (outgoing) particles with spin a in lead j and Aj, Aj are vectors of n real numbers, 
one for each open channel. In terms of incoming (outgoing) creation operator the number 
operators can be expressed as follows 



Eq. dlip can also be recasted in the form 32| : 



XE(X h A~j) = det(I - n E + n E S j A j S A), 



(12) 



(13) 



where I is the unit matrix, is the diagonal matrix of Fermi distribution functions fj{E) for 
particles in the reservoir j and defined as (nE)j a ,jc = fj(E), whereas A is a diagonal matrix 
defined as: (A)^^ = exp(i\j a ). For long measurement times t the total characteristic 
function \ is the product of contributions from different energies, so that 



(14) 



x (\ h X l )=e^ dE kzxB&A) . 

At zero temperature, the statistical average over the Fermi distribution function in Eq. ffTT]) 
simplifies to the expectation value calculated on the state | tp) containing two electrons of 
both spin species for each channel of a given lead up to the energy corresponding to the 
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chemical potential of such lead. Furthermore, in the limit of a small bias voltage V applied 
between the reservoirs, the argument of the integral is energy-independent so that Eq. ([H]) 
can be approximated to 

-. -> r -, -, i M 
Xih, h) - [xo(A T , AjJ (15) 

where only the zero-energy characteristic function appears and M = eVt/h is the average 
number of injected particles. The joint probability distribution for transferring Qi a spin- 
cr electrons in lead 1, Q 2a spin-a electrons in lead 2, etc. is related to the characteristic 
function by the relation (we assume that no polarizers are present): 

P(Qi h Q U , Q2 h -..) = J V dX 1 ^dX n d\ 21 . . . X (h, h) e^-^t e ^Qi . (i 6 ) 

In the rest of the paper we will consider systems where only two counting terminals 
are present. In particular, while the counting terminals are kept at the lowest chemical 
potential, all other terminals are biased at chemical potential eV. For later convenience, we 
write down the most general expression for the characteristic function when spin-cx electrons 
are counted in lead 1 and spin-a' electrons are counted in lead 2: 

Xe(\i*, A 2CT ') = 1 + (e" lA - - 1) (Ng) + (e" tA -' - l) (N^) + 

+ (e~ iAlCT - 1) (e~ iX ^' - 1) (NqN^), (17) 

in the relevant energy range < E < eV. The parameters A corresponding to all others 
terminals are set to zero. 

Using Eqs. (fT5|) . (fT6|) and (IT?]) , at zero temperature, one can calculate the single terminal 
probability distribution: 



M 



M-Qi 



P(Qi„) = ( qI I I | Nq I (18) 

and the joint probability distribution: 

(A/-Qi ct )+(A/-Q 2ct ,) 

P(Qla,Q2a>)= Yl A ™-Q la -Q 2a ,-k B Q la -M+k C k-M+Q 2a , x 

k=Max[M-Q- LtJ ,M-Q 2ly ,} 

x | NgNg 1 | ^) M ~ k f(M, Q la , Q 2a ,,k) (19) 

where A = 1 - (i/j | \ i/j) - (if> \ Ng' \ ip) + (*p\ NgNg 1 \ ip),B = (ip\ N^(l-N^) 
C={^ \ (l-Ng)Ng' \ij>) &ndf{M,Q la ,Q 2a ,,h) = M\/[{k-M + Q 2a ,)\{2M-k-Q la - 
Q 2(7 i)\]. In doing so we have written the expressions for the probability distributions in terms 
of the expectation values of "outgoing" number operators. For Q\ a = Q 2a > = M, Eq. (fT9"|) 
reduces to 

P(Q la = M, Q 2a , = M) = (if, | N^Ng' | V) M • (20) 

When both spin species are counted in one of the terminals the characteristic function is 
different from the one given in Eq. (IT7|) . In particular, the characteristic function for counting 
both spins in terminal 1 reads: 

XsCAi.AjtoO = 1 + (e' lAl - 1) <(A# + (e~ a -' - 1) (iV=f > + 

+ (e" 4Al - 1) {e~ iX ^ - 1) (fe T + N%') + (e" iAl - l) 2 (JV£ T ##) + 

(e~ iAl - l) 2 {e~ iX ^ - l) (NqNqNq u '}. (21) 



where we have set Aij = Aij = Ai. The expression for the joint probability distribution is 
in general complicated, as one can see in Appendix [A] where such expressions for different 
systems are reported. 



III. RESULTS 

The inequality presented in Eq.flSj) can be tested in various multi-terminal mesoscopic 
conductors. In this Section we present several geometries that can be experimentally realized. 
In order to get acquainted with the informations that can be retrieved from Eq.flSj) we start 
from an ideal case in which the entangled pair is generated by some entangler in the same 
spirit as in the works of Refs. @, In Section Till Bl we shall demonstrate that a normal 
beam splitter in the absence of interaction is enough to generate entangled pairs of electrons, 
therefore constituting a simple realization of an entangler. For comparison we also analyze 
the role of superconductivity in creating spin singlets. 



A. Entangled electrons 

In the setup depicted in Fig. [TJ we assume the existence of an entangler that produces 
electron pairs in the Bell state 



a^(E)al(E) ± al(E)al,(E)} | 0) 



"4 J 



3P 



(22) 



of spin triplet (upper sign) or spin singlet (lower sign) in the energy range < E < eV. 
These electrons propagate through the conductors which connect terminals 3 and 4 with 
leads 1 and 2, as though terminals 3 and 4 were kept at a potential eV with respect to 1 and 
2. Our aim is to test the violation of the CH inequality given in Eq. ([5]) for such maximally 
entangled states. 

When the angles 9\ and 82 are parallel to each other, the scattering matrix of the two 
conductors, in the absence of spin mixing processes, can be written as: 



S 



Si 3 
S 2A 



(23) 



where 



S 



13 



r 3 



hi 



( ^3T 



V 



t 



r 3l 





13J 



31T 






\ 

311 



r n J 



(24) 



Here r JfT (tij a ) is the probability amplitude for an incoming particle with spin a from lead j 
to be reflected (transmitted in lead i). For a normal-metallic wire we set = ty^ = vT, 
-VT and rji = Tji = y/l — T, where T is the transmission probability. The 



expression for S24 is written analogously. For simplicity we will assume that S13 and S24 are 
equal. The general scattering matrix relative to non-collinear angles is obtained from S by 
rotating the spin quantization axis independently in the two conductors (note that this is 



S 



possible because the two wires are decoupled). The "rotated" S-matrix is obtained 3l[ by 
the transformation Sg u g 2 = USU\ where U is the rotation matrix given by: 



U 



( u 9l o o\ 

10 
Ue 2 
y I j 



(25) 



where 



cos I 



sin tt 



— sm ~ cos ~ 



(26) 



The probability distributions are now given by the expressions in Eq. (fl8|) and Eq. (fT9|) 
where the state is given by Eq. (|22|) . In the case where both analyzers are present 
we set a — a' =|. The probability distribution when one of the analyzers is removed 
also possesses the structure of Eq. (fl9l since, in this case, the correlators (NqNq) and 
(NqNqNq) appearing in Eq. fpZTT) vanish. In particular when, for example, the upper 
analyzer in Fig. [I] is removed we need to replace N^ 7 with Nq 
For the other correlators one gets: 



7Vq T and Ng' with N$. 



T 



7 1 



and 



T 2 



COS 



^1±^2 



^1±^2 



(27) 

(28) 
(29) 

(30) 



For the single terminal probability distributions in leads i — 1, 2 we get, in the presence 
and in the absence of an analyzer, respectively, 



P 9i {Qi 
P(Qi 



M 
Qi 

M 

Qi 



l - 



M- 



(T) Qi (1 - T) 



M-Qi 



(31) 
(32) 



so that the no-enhancement assumption reads: 

rp\{M-Qi) /-, 



2 



<(1-T) 



(M-Qi) 



1,2 . 



(33) 



Note that the probabilities in Eqs. (1311) and (1321) do not depend on the angles 9\ and 62 
because the expectation values in Eqs. (1271) and (|2"8"|) are invariant under spin rotation. As 
a consequence, the effect of the analyzer is equivalent to a reduction of the transmission 
probability T by a factor of 2, resulting in a shift of the maximum of the distribution. From 
Eq. (l33l) it follows that, for a given number M = eVt/h of entangled pairs generated by the 
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entangler, the no enhancement assumption can be verified only for certain values of T and 
of Qi. This makes clear that the CH inequality of Eq.([Sl) can be tested for violation only 
for appropriate values of M, T and Qi or Q 2 . For example, for a given observation time t 
(i.e. a given M) and a given value of Q, CH inequality can be tested only for transmission 
T less than a maximum value given by the expression 

Qi 



2M-Qi _ i 

T max = — gj- . (34) 

2M- Qi _ 1 

At the edge of the distribution (Qi = M) the no-enhancement assumption is satisfied for 
every T. The window of allowed Qi values where the no-enhancement assumption is satisfied 
gets wider on approaching the tunneling limit. For large M, T max ~ 2 (log The previous 
inequality can be also interpreted as a limit for the allowed measuring time given a setup 
at disposal. Alternatively, given a certain transmission, the no-enhancement assumption is 
verified for points of the distribution such that: 



Qi . log 1 ! 



T/2 



M-log2 + logi^' (35) 

The various probabilities needed to define Sch are collected in Appendix [A] However, it 
is useful to note here that the joint probabilities with a single analyzer are factorized: 

P^(Qi,Q 2 ) = P(Q 1 )P e *(Q 2 ), (36) 

while joint probabilities with two analyzers are not factorized. Furthermore, all such prob- 
abilities have a common factor, T® 1+ ® 2 /2 M , which leads to an exponential suppression for 
large M and Q\ + Q 2 . We shall address the question of whether this also produces a sup- 
pression of Sch in case of violation. 

Let us now analyze the possibility of violation of the CH inequality for different values of 
Qi and Q 2 . First consider the situation where the entangler emits a single entangled pair of 
electrons in which case P fll '" 2 (l, 1) = | N$ N$ | tp), P~^ 2 (l, 1) = (ip\ (n£+Nq 1 )N% | ip) 
and P dl ~(l, 1) = (if) ] N$(N$ + N~§) \ tp). We find that the CH inequality is maximally 
violated for the following choice of angles: 6 2 — Q\ = 8' 2 — 9[ = 3n/4. More precisely we 
obtain: 

Sch = T 2 —^— (37) 



which is equal to the result obtain for an entangled pair of photons [30(, where T plays 
the role of the quantum efficiency of the photon detectors. In the more general case of 
Q 1 = Q 2 = M, for M > 1, we have 



P dl ' d2 (M, M) 



2 



2M r 2 fe 1 ±e 2 ^ lM 

sin 



M 



2 



rp 2 M 

P e ^-(M,M) = p-> d2 (M,M) = (38) 
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so that the no-enhancement assumption is always satisfied and the quantity Sch can be 
easily evaluated: 



T 



2M 



CH 



2 M 



sin 2M - sin 2M ^ + sin ^ °±±h + ^2 



(39) 



The rotational invariance makes P dl ' and P ,02 independent of angles, and P 01 ' 02 dependent 
on the angles through 1 ^ 2 . This allows us, without loss of generality, to define an angle 
such that 29 = 9 1 ± 9 2 = 0[ ±9 2 = 9[ ±9' 2 = (#1 ± 9' 2 ) /3. As a result Eq.© takes the form: 



<S OT = 3P^ 2 (Q ll Q 2 ) - Pjf (Qi, Q 2 ) - Pi,-(Qi, Q 2 ) - ^, 2 (Qi, Q 2 ) < (40) 

where P® 2 = P 6 * 1 ' 6 * 2 and Pi_ = P 01 ~. It is useful to define the reduced quantity Sch = 
Sqh / {T 2M /2 M ) which is plotted in Fig. as a function of G for different values of M 
(note that since P 6l '"(M,M) = (T 2M /2 M ), S CH is nothing but S C h/P 9u ~(M, M)). The 
violation occurs for every value of M in a range of angles around = 7r/2 (note that <Sch 
is symmetric with respect to vr/2). The range of angles for which Sch is positive shrinks 
with increasing M, while the maximum value of Sch decreases very weakly with M (more 
precisely, S C h ^ This means that the effect of the factor T 2M /2 M on the value 

of Sch is exponentially strong, making the violation of the CH inequality exponentially 
difficult to detect for large M and Qi = Q 2 = M. The weakening of the violation is mainly 
due to the suppression of the joint probabilities. As we shall show later, by optimizing all 
the parameters it is yet possible to eliminate this exponential suppression. 

Let us now consider the violation of the CH inequality as a function of the transmitted 
charges. We notice that the CH inequality is not violated for the off-diagonal terms of the 
distributions (when Q\ 7^ Q2), meaning that one really needs to look at "coincidences". 
Therefore we discuss the case Qi = Q 2 = Q < M (remember that the no-enhancement 
assumption is satisfied only for T < T max (Q)). In Fig. [3] we plot the quantity Sch for 
M = 20 as a function of and different values of Q. The transmission T is fixed at the 
highest allowed value by the no-enhancement assumption, which corresponds to the smallest 
Q considered T max (Q = 1) = 0.06917. Fig. [3] shows that the largest positive value of Sch 
and the widest range of angles corresponding to positive Sch occur for Q = 1, i.e. for a 
joint probability relative to the detection of a single pair. One should not conclude that, in 
order to detect the violation of the CH inequality, only very small values of the transmitted 
charge should be taken. We have in fact considered T = T max relative to Q = 1 and the 
maximum violation, for given M and Q, always occurs at T = T max - In order to get the 
largest violation of the CH inequality at a given M and Q one could, in principle, choose 
the highest allowed value of T for each value of Q (T — T max (Q)). We show in Fig. |4]the 
corresponding plot, to be compared with Fig. [3J For every Q < M the violation occurs in the 
same range of angles, namely 7r/4 < < ir/2, because of the following properties of the joint 
probability distributions: Pf 2 {QuQ 2 ) = Pf%{Q l ,Q 2 ) = P 1 _(Q 1 ,Q 2 ) for = tt/4. This 

implies that S CH (Q = tt/4) = 0, and P^(Qx',Q 2 ) > P${Qu Q2), P£-{Qu Q*)> P-AQi, Q2) 
for 7r/4 < < it/2. Furthermore, in this specific case of M = 20, we find that the maximum 
values of S occurs at Q = 8. 

In Fig. Owe plot the maximum value of S, with respect to and T, as a function of Q for 
different values of M. Several observations are in order. For increasing M, the position of 
the maximum, Qmax is very weakly dependent on M. Remarkably, the value of the maximum 
of the curves does not decreases exponentially, but rather as 1/M 2 . Despite the exponential 
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suppression of the joint probability with M, the extent of the maximal violation scales with 
M much slowly (polynomially). 

It may be useful to look at the same situation from a different perspective. Given a 
certain transmission T (i.e. fixing the transport properties of the conductors) we want to 
find when the CH inequality is maximally violated. For a given observation time t, the 
no-enhancement assumption Eq. (1341) imposes a minimum value for Q. In Fig. we plot 
the quantity Sch-, maximized over the angle G and Q, as a function of T for different M. 
The curves are piecewise increasing function of T, where the discontinuities correspond to 
an increase of the value of Q by one imposed by the no-enhancement assumption. More 
precisely, when T is increased above a threshold for which Eq. (]34p is not satisfied, one 
needs to increase Q by one unit in order for this condition to be recovered. The result of 
this is a jump in the values of the probabilities that leads to a discontinuity of Sqh- Fig. [6] 
allows to choose the best values of M and Q to get the maximum violation. 

If the entangler is substituted with a source that emits factorized states, the CH inequality 
given in Eq.([3J is never violated. In this case, in contrast to Eq. (1221) . the state emitted by 
the source reads: | ip) = aLaL | 0) . All the previous calculations can be repeated and 
we find, as expected, that the characteristic functions factorizes, so that the two terminal 
joint probability distributions are given by the product of the single terminal probability 
distributions. 



B. Normal beam splitter 

We are now ready to analyze realistic structures by replacing the shaded block in Fig. [1] 
(which represents the entangler) with a certain system, and discuss the CH inequality along 
the lines of Section IHI Al We first consider a normal beam splitter (shaded block in Fig. [7]) 
in which lead 3 is kept at a potential eV and leads 1 and 2 are grounded so that the same bias 
voltage is established between 3 and 1, and 3 and 2. The two conductors, which connect 
the beam splitter to the leads 1 and 2, are assumed to be normal-metallic and perfectly 
transmissive, so that the S-matrix of the system for 9\ = 62 = is equal to the S-matrix of 



the beam splitter, which reads [40 



-(a + b) y/e y/e 
S = I y/l a b I . (41) 
y/e b a 



In this parametrization of a symmetric beam splitter a = ±(1 + yl — 2e)/2, b = =(=(1 — 
\/l — 2e)/2 and < e < 1/2. For arbitrary angles 9\ and 82, the S-matrix is obtained rotating 
the quantization axis in the two conductors independently by applying the transformation 
Se 1 fi 2 = USW , where hi is the rotation matrix given by: 



I 






















Ue 2 



U = U 6l (42) 
\0 U e J 

and Uq is defined in Eq. ( 1261) . This procedure is valid as long as no back scattering is present 
in the conductors. The probability distributions are given by Eqs. ffl8l) and ffl9l) where the 
state I ip) is now factorisable: 

|^) = a t 1T (E)a t u (E) |0> (43) 
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in the energy range < E < eV. Analogously to what was done in Section IIII Al when 
both analyzers are present we set o = a' =f. When only one analyzer is present, however, 
one has to use the correct characteristic function of Eq. (l2Tj) . since one of the two additional 
correlators does not vanish. Namely, {Nq Nq) = e 2 and {N Q N Q Nq) = 0, when the upper 
analyzer, for example, in Fig. [7] is removed. For the other expectation values we get: 

^\N^ \^) = (^\NS \^) = e, (44) 

ty\N#\1>)=e, (45) 

(^KiV?^)= e 2 sm 2 (^) (46) 

and 

M^M=eW(^Y (47) 



obtaining the joint probability distributions reported in Appendix |A] The above number 
operator expectation values are equal to the case of the entangler when e is replaced by 
T/2, whereas the cross-terminal correlators are equal in the two cases if e is replaced with 
T/y/2. From this follows that the characteristic functions for the beam splitter possess the 
same dependence on the angle difference as the corresponding characteristic functions for 
the entangler (Section IIII Al) but have a different structure as far as scattering probabilities 
are concerned. In particular, as expected [32J, the cross-correlations vanish when the two 
angles are equal. On the contrary, when the angle difference is ir cross-correlations are 
maximized. Furthermore, when only one analyzer is present the characteristic function 
shows no dependence on the angle, but it is not factorisable, in contrast to the case of the 
entangler. As a result, the single terminal probabilities, given by Eq. (fT8]) . are equal in the 
two cases provided that e is replaced with T/2. The joint probabilities for Q\ = Q2 = M 
are equal in the two cases if e is replaced with T j\[2 (however, this replacement is not valid 
in general for joint probabilities with Qi,Q 2 7^ M): 



701,02 



(M, M) 



2 ■ 2 
e sin 



(48) 



P 6l '-(M,M) = e 2M (49) 
The no-enhancement assumption is verified when 

Q 

1 2m=q _l 

6 < -— -, (50) 

^ 2^^q — I 

which equals the condition of Eq. (l34"j) once e is replaced with T/2. Let us first consider 
the case for which Q\ = Q 2 = M. We obtain an important result: the CH inequality is 
violated for the same set of angles found for the case of the entangler, although to a lesser 
extent, since the prefactors in Eqs. (1481) and (1491) now varies in the range < e 2M < 
In particular, in the simplest case of M = 1, corresponding to injecting a single pair of 
electrons, the maximum violation corresponds to Sch = ^ 4 ~ 1 , which is a half of the value 
for the entangler. Furthermore, the plot in Fig. [2] is also valid in the present case with Sch 
defined as Sch = Scn/ e2M , by replacing with e. This means that a geometry 
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like that of the beam splitter enables to detect violation of CH inequality without any need 
to resort to interaction processes to produce entanglement. 

Also here we consider the case for which Qi = Q2 = Q < M, where interesting differences 
with respect to the case of the entangler are found, i) We find that the violation of the CH 
inequality is in general weaker, meaning that the absolute maximum value of Sqh is smaller 
than in the ideal case of the entangler. ii) The weakening of the violation with increasing M 
is determined by the suppression of the probability by the prefactor (g^Qi+Q 2 , Remarkably, 
the maximum value of 5 max decreases like 1/M, therefore even slower than for the ideal case, 
iii) Violations occur only for values of Q close to 1, even for large values of M: to search 
for violations one has to look at single- or few-pair probabilities and therefore, because of 
the no-enhancement assumption, to small transmissions e. iv) Interestingly, for Q = 1 the 
quantity Sch is positive for any angles, although the largest values correspond to close 
to 7r/2 (see Fig. [H]). We do not find any relevant variation, with respect to the discussion in 
paragraph IIII Al for probabilities relative to Qi 7^ Q 2 . 

It is easy to convince oneself that the final state calculated from the initial one f)43p using 
the S-matrix (]41[) contains an entangled part. In Ref. [4l[ this fact was already noticed, but 
for an incident state composed by a single pair of particles impinging from the two entering 
arms of a beam splitter. For mesoscopic conductors, entanglement without interaction for 
electrons injected from a Fermi sea has been discussed by Beenakker et al 19l . 



C. Superconducting beam splitter 

In many proposals superconductivity has been identified as a key ingredient for the cre- 
ation of entangled pairs of electrons. The idea is to extract the two electrons which compose 
a Cooper pair (a pair of spin-entangled electrons) from two spatially separated terminals. 
Here we showed that it is not necessary to have superconducting correlations. Nevertheless, 
in view of the recent interest in entanglement created by pairing correlations, it is useful to 
analyze also the case of a superconducting beam splitter 0, HI depicted in Fig. which 



consists of a superconducting lead (with condensate chemical potential equal to /1) in contact 
with two normal wires. The wires are then connected to two leads attached to reservoirs 
kept at zero potential. This is basically what is obtained by replacing the entangler of Fig. [1] 
by a superconducting lead with two terminals. 

The system can be decomposed into two subsystems: on the left-hand-side of Fig. [9] we 
place the superconducting slab attached two normal terminals (5 and 6) characterized by a 
reflection amplitudes matrix R' s defined, in terms of the particle operators, by: 

4> jaa (E) = E E E K(£)W,^ *wAE). (51) 

k=5,6 P=e,h<r'=1,l 

Here j — 5,6 and the additional indexes a and (3 refer to the particle-hole degree of freedom, 
in particular a = e for particles and a = h for holes and [. . .]j aa ,kf3a' represents the specified 
element of the matrix. Note that R' s is block diagonal in spin indexes so that 
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with 



5eT 
5hl 



&5hl 

a 6e| 
V ®6H ) 



( Pee 
Phe 



Pph 
Phh 



T' T, 



ell 



\ T he T hh 



Tee 
The 
Ppp 
Php 



Teh \ 
Thh 
Peh 
Phh ) 



(53) 



where p ee (phh) is the normal reflection amplitude for particles (holes) in terminal 5, p e h 
(phe) is the Andreev reflection for a hole to evolve into a particle (particle to evolve into a 
hole) in terminal 5. r ee (Thh) is the normal transmission amplitude for particles (holes) to be 
transmitted from terminal 5 to terminal 6, r e h (The) is the Andreev transmission amplitude 
for holes (particles) in terminal 5 to be transmitted in terminal 6 as particles (holes). Primed 
amplitudes refer to reflections occurring in lead 6 and transmissions from lead 6 to lead 5. 

On the right-hand-side of Fig. [9] we have the subsystem composed of two identical decou- 
pled conductors characterized by the 16 x 16 scattering matrix 



1 R c T c 
' 1 T r R' 



(54) 



The four submatrices in Eq. flMl) are block diagonal in spin space, for example R c can be 
written as: 

Rl 



Rc 



where R\ is a diagonal matrix defined by 



Rl 



(55) 



/ 03eT \ 



Rl 



\ <pAhl J 



( fl3eT \ 
®3hl 
a>4ef 
\ 0>4hl J 



Rl 



/r 3eT \ 

r 3H 

r 4eT 

r 4H ) 



(56) 



R\ is defined like R\ exchanging j with j, whereas T° is defined similarly to R° replacing 
r 3aa with ti aa and r iaa . with t 2a(T . The matrices R'° and T C /<J are defined analogously using 
the amplitudes r laa , r 2aa , t' laa and t' 2aa . The spin quantization axis of the two wires can be 
rotated independently as in paragraph IIII Al by applying the transformation Sg 1 1 $ 2 = US C U\ 
where U is defined in Eq. ( |25|) . obtaining the scattering matrix 



5 ei,0 2 



T c Rr 



(57) 



The overall matrix of reflection amplitudes is calculated by composing the scattering matrices 
relative to the two subsystems |42|: 



Rf nt Rc ~l~ T c 



I - R' a R c 



r't' 



where R' tQt is defined by 



aa,k/3r 



afc/3r(-E , ) ! 



(58) 



(59) 



fc=l,2/3=e,/iT=U 
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with j running from 1 to 2. The characteristic function can now be calculated through 
Eq. (fT3|) taking R' tot (E) as scattering matrix. In the present case, where superconductivity 
is present, the diagonal matrix of Fermi distribution functions is defined as [nE]jacr,jacr = 
fja(E), fja{E) = [1 + exp( ; |^)]- 1 and [A] j(Wtjaa = exp(ia\ ja ), with j = 1,2. By choosing 
Ai| = A21 = we achieve the goal of counting excitations with spin- up component. The 
case where one of the analyzers is removed, for example in lead 1, is implemented by setting 
Ai| = Ai| = Ai and 61 = 0, i.e. by counting electrons in lead 1 regardless their spin. 

In the limit of zero temperature and small bias voltage, we only need the scattering 
amplitudes at the zero energy (Fermi level) so that the overall characteristic function can 
be approximated like in Eq.( jT5l) . We parametrize the matrix S c of the wires as follows: 

r-3ea = r 4 ea = \A ~ T, T\ e(J = T 2 ea = V 1 = T, t lea = t lea = Vf and t' lea = t' 2ea = ~VT, 

where T is the wire transmission probability of the wires. The amplitudes relative to hole 
degree of freedom are determined from the ones above by making use of the particle-hole 
symmetry. 

Although Andreev processes are fundamental for the injection of Cooper pairs, in the case 
where Andreev transmissions only are non-zero and T = 1 the joint probabilities factorize 
in a trivial way 

P u 2 (Qi> Q2) — 8qi,2m8q 2 ,2m P dl ' {Qi, Q2) — <5<2i,2m£q 2 ,4M , (60) 

in such a way that the CH inequality is never violated. This apparent contradiction is due 
to the fact that in this situation the scattering processes occur with unit probability, so 
that the condition of locality is fulfilled. Non- locality can be achieved by imposing T < 1. 
In the limit T C 1 we obtain the probabilities P 1, 2 {Q\, Q2) and P~ ,e2 (Qi, Q2) reported, 
respectively, in Eqs. (AED and (AED of the Appendix, which reduce to 

P e ^(M,M) = 

and 

M 

(62) 

for Q 2 = Q 3 = M, with A = 1 + T he r£ e . Eqs. (EE)) and ([62]) are equal to Eqs. (JHEj), relative 
to the case of an entangler, once 2T 2 A 6 /[A — T(A — l)] 8 is replaced with T 2 /2. From this 
follows that superconductivity leads to violation of the CH inequality. For A = 2, i.e. perfect 
Andreev transmission, the quantity 2T 2 A 6 /[A — T(A — l)] 8 tends to T 2 /2 in the limit T — ► 
so that the analysis of Section Till Bl relative to the case Q\ = Q2 = M applies also here. 

IV. CONCLUSIONS 

In mesoscopic multiterminal conductors it is possible to observe violations of locality 
in the whole distribution of the transmitted electrons. In this paper we have derived and 
discussed the CH inequality for the full counting electron statistics. In an idealized situation 
in which one supposes the existence of an entangler, we have found that the CH inequality 
is violated for joint probabilities relative to an equal number of electrons that have passed 
in different terminals. This is related to the intuition that any violation is lost in absence 
of coincidence measurements. The extent of the violation is suppressed for increasing M 



2T 2 A« 



[A — T(A — 1)Y 



IVl 



sm 



9i + 9 2 



ivi 



(61) 



,02 



(M, M) 



2T 2 A* 
[A - T(A - 1W 
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(average number of injected pairs), however such a suppression does not scale exponentially 
with M like the probability, but instead decreases like 1/M 2 . This means that the detection 
of violation does not become exponentially difficult with increasing M. For fixed transport 
properties we analyzed the conditions, in terms of M and number of counted electrons, for 
maximizing the violation of the CH inequality. 

The violation of the CH inequality could be achieved in an experiment. Indeed we tested 
the CH inequality for two different realistic systems, namely a normal beam splitter and 
a superconducting beam splitter. Interestingly we find a violation even for the normal 
system, even though weaker with respect to the idealized case of the entangler. In this case 
the violation is again suppressed for increasing observation time, but scales like 1/M. We 
analyzed the superconducting case in the limit of small transmissivity and we also find a 
violation of the CH inequality to the same extent with respect to the case of the entangler. 

It is important to notice that the analyzers should not affect the scattering properties 
of the system as in the case of ferromagnetic electrodes. In the latter case, in fact, the 
probability density of the local hidden variables would also depend on the angles 81 and 9 2 . 

We believe that the results derived in this work may be of interest for the understanding 
of the statistics of electrons in mesoscopic conductors. It is however important to look for 
experimental tests of our claims. In this respect two possible schemes for measuring the 
counting statistics have been recently proposed in Ref. [43j ]. Since solid state devices are 
considered promising implementations for quantum computational protocols, this line of 
research does not seem interesting only from a fundamental point of view, but may be of 
clear relevance for the actual realization of solid state computers. 
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Appendix A: PROBABILITY DISTRIBUTIONS 

In this appendix we give the general expressions for the joint probability distributions 
used in the paper to discuss the CH inequality. 

1. Entangler 

In the case of an entangler we find 



P ei >-(Qi,Q2) 



P^ 2 (Qi,Q 2 ) 




M-Q 2 



M-Q 2 



(lb) 



(la) 
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and 



Min[Qi+Q 2 ,M] 

E 

fc=Max[Qi,Q 2 ] 

rp{Ql+Q2) 



M 



k 



k J \2k - Qx - g 2 



2k-Q 1 -Q 2 



2 M 

1 — T sin 



2(1 — i ) + 1 sin 



/ 01 ± 2 



i 2fc-Qi- 



sm 



k-Q 2 

~l M-fc 
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^1 ±^2 



(2) 



2. Normal beam splitter 

The joint probability P ,62 (Qi,Q 2 ) used in Section Ull Bl is 

Mm[(M-Q 1 )+(M-Q 2 ),M] 

P ei ' e %Qi,Q2) = 



x e 



x 



fc=Max[M-Qi,M-Q 2 ] 
(Q1+Q2) 



M 



A- 
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/ 01 — 00 
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sm 



01-02 
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The single-analyzer joint probability P ,e2 {Qi,Q2) reads: 

Qi Min[A/-fc,Q 2 ] 

E E 

fc=0 «=Max[0,(Qi-fc)+(Q 2 -fc)] 

'k + l-Q 2 
Qi-k 



M\fM-k\f k 
k)l I )\k + l-Q 2 



(3) 



x 



(* q[ - k 2 ) ~ 3e + 2e2 ] M_fc "' [1 ~ e]l[2 ~ 3e] 



2k+l-Q 1 -Q 2 



(4) 



with < Qi < 2M and < Q 2 < M (note that the sum on I has to be performed only 
when the lower limit is less than or equal to the upper limit). 



3. Superconducting beam splitter 

The joint probability P u °*(Qi, Q 2 ) used in Section Ull CI is 

Um[Q 1 +Q 2 ,M] 

P ei ' 62 (Qx,Q2) = £ 



fc=Max[Q!,C 
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where A = 1 + r hp r^ p . 

The single- analyzer joint probability P~ ,62 (Qi-,Q2) reads: 

^(ft.fc) = (*) ( „_ r g_ 1)18 ) M [1 -4T + 6Tr- (6) 

for Q 1 > Q 2 and P-^{Qx,Q 2 ) = for Q 1 < Q 2 . 
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Figure 1: Idealized setup for testing the CH inequality for electrons in a solid state environment. It 
consists of two parts: an entangler (shaded block) that produces pairs of spin entangled electrons 
exiting from terminals 3 and 4. These terminals are connected to leads 1 and 2 through two 
conductors described by scattering matrices S13 and S24. Electron counting is performed in leads 
1 and 2 along the local spin-quantization axis oriented at angles Q\ and 62- 
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Figure 2: The quantity S C h = S C h/{T 2M /2 m ) is plotted M S iX function of the angle for different 
numbers M of injected entangled pairs by the entangler. The range of angles relative to positive 
values shrinks with increasing M, while the value of the maximum slightly decreases. 
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Figure 3: The quantity Sch is plotted as a function of the angle G for M = 20 and T = 0.06917, 
which corresponds to the highest value allowed by the no-enhancement assumption for Q = 1. The 
curves are relative to different values of Q = [1,4]. Note that for Q > 4 the variation of Sch over 
the whole range of is small on the scale of the plot. Violations are found only for Q = 1 and 
Q = 20. 
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Figure 4: The quantity Sch is plotted as a function of the angle for M = 20 and T set to 
the highest value allowed by the no-enhancement assumption, different from each Q. The curves 
are relative to different values of Q = [1,20]. The maximum of Sch increases with Q reaching 
its largest value for Q = 8 and decreasing for Q > 8. Note that the variation of Sch with for 
Q = 20 it is not appreciable on this scale. 
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Figure 5: The maximum value of the quantity Sqh, evaluated over angles G and transmission 
probabilities T, is plotted as a function of Q. The curves are relative to different values of M 
ranging from 10 to 30. For points corresponding to the maximum of the curves we indicate the 
corresponding value of transmission T. 
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Figure 6: The maximum value of the quantity Sqh, evaluated over angles and number of counted 
electrons Q, is plotted as a function of T. Both curves, relative to M = 10 and M = 20, exhibit 
discontinuities which correspond to an increase of the value of Q by one. This increase is imposed 
by the no-enhancement assumption, Eq. ([35]) . which depends on the value of T. We indicate the 
value of Q which corresponds to the largest violation. 
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Figure 7: Setup of a realistic system consisting of a normal beam splitter (shaded region) for testing 
the CH inequality. Bold lines represent two conductors of unit transmission probability. A bias 
voltage equal to eV is set between terminals 3 and 1 and terminals 3 and 2. 
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Figure 8: The quantity Sqh for a normal beam splitter is plotted as a function of the angle for 
three values of M = eVt/h = 10, 20, 100 when Q = 1. Interestingly, Sqh is positive for every angle 
and its maximum value decreases like 1/M. 
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Figure 9: Setup of a realistic system consisting of a superconducting beam splitter (shaded region) 
for testing the CH inequality. Bold lines represent two conductors of transmission probability T. 
The superconducting condensate electrochemical potential is set to /j,, while terminals 1 and 2 are 
grounded. 
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